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Methods of phase measurements at VLF are briefly described. Some results of diurnal 
phase measurements, sunrise effects, and solar flare effects at medium distances and northern 
geographical latitudes are summarized. Two inhomogeneous and anisotropic models are 
used for the interpretation of the measurements: the first corresponding to the daytime 
ionosphere and the second corresponding to the nighttime ionosphere. The models are 
related to homogeneous models with sharp boundaries and effective parameters of conductiv- 
ity, dielectric constant, and equivalent height. The diurnal phase lag, depending on dis- 
tance, and the sunrise effects at 1000 km and 2000 km are interpreted by these models. 
Two critical distances exist where the nighttime phase gains a phase advantage of 360 degrees 
relative to the daytime phase. An electron profile of the daytime lower ionosphere is deduced 
from measured daytime variations of equivalent height. The strength of ionizing solar 
radiation during a solar flare and the effective attachment coefficient of the lower ionosphere 
are derived from measured variations of the equivalent height during a solar flare. 



1. Introduction 

VLF waves (3 to 30 kc/s) are increasingly used as 
communication medium for lime and frequency 
comparison and for long range navigation systems. 
In all these cases, the frequency of a cw-transmitter 
is compared with the frequency of an oscillator at 
the recieving station, which is independent of the 
transmitter oscillator. 

VLF wave propagation fakes place in the wave- 
guide between earth's surface and the lower iono- 
sphere. This waveguide behaves as a transformer 
and depends on frequency, time, and location. Each 
change of the transformer characteristics of the wave 
guide is accompanied by a change of phase of the 
transmitter at the receiving station relative to the 
receiver oscillator. Therefore, the transformer char- 
acteristics and their variations must be known in any 
practical application. 

In the atmospheric waveguide primary, the re- 
flection characteristics and the equivalent height of 
the lower ionosphere vary. On the other hand, the 
influence of the inhomogeneous earth's surface, with 
different coverings in winter and summer, is of subor- 
dinate importance. Therefore, VLF waves are a 
very useful tool for studying the lower ionosphere. 

Observations of VLF propagation at medium dis- 
tances are especially suitable for this research because 
the strong influence of the skywaves increases as to 
distance increases while the earth's magnetic field 
has only a small effect. Otherwise, the region of 
propagation behaves nearly as a space which is 
homogeneous in the horizontal direction. There- 
fore, it is possible to use a relatively simple model 
for the theoretical interpretation. 



Medium distances shall be distances where the 
propagation of electromagnetic energy does not take 
place predominantly by one mode. At daytime 
this range comprises of distances between about 300 
and 3000 km. Within 300 km the principal wave 
is the groundwave. Beyond 3000 km the first wave- 
guide mode is the chief carrier of energy. At night- 
time the range becomes larger because of the smaller 
absorption of the ionospheric skywaves. 

A great deal of early work of VLF research lias 
been done by the Cambridge group. The results 
are summarized by Bracewell et at., [1951]. The 
easiest method of phase measurement is to lock the 
receiver oscillator with the transmitter oscillator, 
e.g., by telephone circuits [Bracewell, 1952]. Natu- 
rally this is possibly only for relatively short dis- 
tances and suitable local conditions. Pierce [1955, 
1957] compared the frequency of an independent 
oscillator in Cambridge (Mass.) with the highly 
stabilized frequency of GBR (Rugby), 16 kc/s, and 
first realized useful phase measurements over large 
distances. 

This paper reports results of diurnal phase meas- 
urements at medium distances and their interpreta- 
tion. Attention is restricted to regular changes 
including sunrise effects and solar flare effects, but 
does not discuss nighttime anomalies like magnetic 
storm effects or cosmic ray effects. The reason, 
therefore, is that the nighttime ionosphere is very 
unstable and becomes highly inhomogeneous in the 
horizontal direction at such anomalies, consequently 
a quantitative interpretation is extremely difficult. 
Since only the phase together with the amplitude of 
the field strength gives unique information, it is 
sometimes inevitable for the interpretation of the 
measurements to use the amplitude too. 
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2. Results of Phase Measurements 

2.1. Method of Phase Measurement 

The phase of a transmitter at the receiving station 
has the form 



ti=wt — kp-\-<p(p,t) 



(i) 



where <p(p, t) is a variable form which depends on both 
range, p, and time, t. The receiver oscillator has 
the phase 

Q E =u E t J r<p E . (2) 

Here o> and u E are the angular frequencies of the 
transmitter and the receiver oscillator, respectively; 
t is the time; k=w/c is the wave number; c is the 
velocity of light; and p is the distance between 
transmitter and receiver. The phase 



cot — kp=u(t—p/c) 



(3) 



is the phase of an undisturbed wave propagating in 
free space; <p(p, t) is the phase difference from this 
ideal value; and <p E is an arbitrary phase constant of 
the receiver oscillator. 

If the two frequencies are mixed in a converter, a 
difference frequency will be obtained, and the result- 
ing phase difference is 



12 — 12#=(co — oo E )t-\-<p(p, t) — kp — <p E . 



(4) 



Since the oscillators of the transmitter and the re- 
ceiver generally are not synchronized, the result of 
phase measurement is undecided by the factor 



kp+(p E . 



(5) 



Furthermore, each frequency difference of the two 
oscillators gives a linear phase change with time 
which adds to the diurnal variation <p(p, t). Often 
the frequency of the receiver oscillator changes by 



o) E =o) EQ J rat. 



(6) 



Then the phase change with time even is of parabola 
form. Such an oscillator run can be compensated 
[Utech and Zahn, 1961]. 

The maximal diurnal change of <p is of the order 



kday - knight |—27T. 



(7) 



In order to measure primary, the change of <p it 
must yield within one daily interval (T^IO 5 sec) 



or 



{o>-oo e )T«2t 

A/«_L 



At the frequency /= 10 kc/s we get 

4L 



/ 



«io- 9 . 



(8) 



(9) 



At=— > 



(11) 



is used when measurements at different frequencies 
are compared. Finally, the phase lag is defined by 

—A<p=<p midd&y —(p. (12) 



2.2. Results of Diurnal Phase Measurements 

Figure 1 shows results of phase measurements of 
GBR (16 kc/s) in Berlin for the year 1962. The 
distance between transmitter and receiver is 980 km. 
The monthly means of hourly median values of A<p 
are scaled versus time of day. The propagation path 
is directed nearly from west to east. The shift of 
the daytime maximum and the shift of the night- 
time minimum to earlier hours than midday and 
midnight, respectively, are striking. The difference 
time of sunrise or sunset at the locations of the 
transmitter and receiver is smaller than 1 hr. The 
variation of phase must be caused, therefore, by a 
smooth variation of the equivalent reflection height 
and the reflection characteristics of the lower iono- 
sphere. On the other hand, Pierce [1957] and Chilton 
et al. [1962], find a trapezium shape of the diurnal 
phase variation for greater distances. The duration 
of the slope of the trapezium depends on the time 
during which parts of the propagation path lie in 
twilight. The duration becomes great for propaga- 
tion paths directed from east to west or vice versa 
and are small for paths directed from south to north 
or vice versa. 

If only the first mode is the carrier of energy (for 
distances greater than about 3000 km), it is sufficient 
for an approximate interpretation of the transition 
to assume that the phase of the mode is propagating 
with a constant phase velocity at daytime and a 
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The frequencies of the oscillator must be stabilized 
at least for 1 in 10 10 , if a useful diurnal phase measure- 
ment at VLF shall be possible. The method of 
achieving this extreme condition is beyond the scope 
of this paper. Details may be found in Pierce [1957], 
Crombie et al. [1958], and Eppen and Heydt [1959]. 

An accuracy of 1 in 10 10 of the oscillator frequency 
is already the upper limit which can be realized 
within economic limitations. The oscillator of GBR 
(16 kc/s) has this accuracy [Pierce et al., I960]. 
Therefore, the determination of the reference phase 
is not unique. Yet, the experiments show that the 
lower ionosphere at daytime is rather stable, and its 
behavior at midday only is changing a little from day 
to day. This is the reason that the phase value at 
midday will be used as reference phase. 

The following data will show the phase difference 

A(f=(f— ^midday, (10) 

although for some purposes the time delay, 




12 74 76 18 20 22 24 02 04 06 08 10 12 

— time [GMT] 

Figure 1. Monthly means of diurnal -phase difference A<p 
of GBR versus time observed in Berlin during the year 

1962. 

The dashed lines indicate time of sunrise (SR) and time of sunset (SS), the dash- 
dotted line indicates time of midnight at the middle of the propagation 
path. 

different constant phase velocity at nighttime. In 
regions of transition from night to day, the portions 
of phase paths are simply added. For example, if p 
is the percentage portion of the nighttime phase path, 
t hen the whole phase at the receiver is [Crombie et al., 
1958] 

"knight), (13) 



(14) 
(15) 



is the measured maximal phase difference relative 
to the phase at midday. Since the change of p 
is approximately linear, (14) gives a trapezium 
variation. 



il=o)t- 


-(1— p)(kp— ^day)— V 


and thus 






&<p=pA<p mj 


where 






&<Pm — knight <Pday 




distance g [km] 

Figure 2. Calculated maximal diurnal phase lag — A^> m 
versus distance of a 16 kc/s transmitter (full lines). 

For details see text. The circles are measured 16 kc/s values, the crosses are 
measured 18 kc/s values. 

With help of (14), Blackband [1962] and Hampton 
and Hill [1962] were able to explain the daily phase 
variations over propagation paths crossing arctic 
regions. In winter or summer there exist two ex- 
tremes of daily phase variation, because parts of the 
propagation paths are dark and illuminated, respec- 
tively, during the whole day. For shorter distances, 
the propagation must be interpreted by ray theory. 
Here it is better to describe the transition from night 
to day by a steady change of the ionospheric char- 
acteristics over the whole propagation path (see sec. 
3.4). 

Casselman et al., [1959], Blackband [1962b], and 
Westfall [1961] have reported the maximum phase 
lag — Acp m for different distances from the transmitter 
and different frequencies between 10 and 20 kc/s. 
Some of these results for the frequencies 16 and 18 
kc/s together with the data from figure 1 and 3 are 
plotted versus distance p in figure 2. 

Contrary to the presentation of these data by 
Blackband, no attempt has been made to eliminate 
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Figure 3. Diurnal phase difference Ap of GBR versus time measured in 

Lindau (Germany). 



the groundwave at smaller distances. For 18 kc/s, 
only the right scale — Ar m is significant. In all cases 
the vertical component of the field strength E z has 
been measured. The strong minimum of phase lag 
in the region between 2000 and 2500 km is striking. 

For all distances, the measured values Ar m are 
negative. A remarkable exception from this rule has 
been observed in Lindau (Germany) . Figure 3 shows 
two successive daily phase variations of GBR. On 
the first day there exists the normal behavior of vari- 
ation (plotted in fig. 2) which fits to the other values. 
The next day an inversion of phase lag is observed. 
The maximum difference of phase on both days is 
about 360°. Such phase inversions are rare and 
have been observed up to now only at distances 
round 700 km from the transmitter. 

Blackband [1962b] has shown that for great dis- 
tances the value — Ar w /p is nearly constant at con- 
stant frequency and decreases with frequency down 
to a minimum at about 16 kc/s. In figure 4 the 
measured values are plotted versus frequency. 
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Figure 4. Calculated values — Ar m /p versus frequency (dashed 

line). 



The points are measured values according to Blackband [1962b]. 



Blackband did not find any significant dependence 
of this data on direction of propagation (influence of 
earth's magnetic field) or propagation path (influence 
of the distribution of continents and oceans). 

2.3. Sunrise Effects 

At medium distances the sunrise effects are very 
regular and marked in phase and amplitude of VLF. 
The sunset effects are much weaker and not so regu- 
lar. As a result of the averaging in figure 1, the 
sunrise effects are smoothed out. Two examples of 
sunrise effects of GBR, measured in Berlin, are 
shown in figure 5. The ordinates are the monthly 
means of July and December 1962, respectively, 
with phase and amplitude separated. The abscissa 
is the time difference centered to the time of zenith 
distance x = 96° at the middle of the propagation 
path. That is equivalent with the time of sunrise 
in 35 km height above the middle of the path. In 
summer a deep minimum of amplitude and a turning 
point of phase can be seen; in winter there is a maxi- 
mum followed by a minimum in amplitude and like- 
wise a turning point in phase. The beginning of the 
sunrise effect is at about x = 99°. This corresponds 
to sunrise at the height of 80 km above the middle of 
the path. At 2000 km the phase of sunrise effect 
can surpass the daily variation as Blackband [1961] 
has shown (see fig. 13). 

2.4. Solar Flare Effects 

A solar flare lowers the whole daytime portion of 
the lower ionosphere. The homogeneous conditions, 
therefore, will not change during a solar flare, and 
the same theoretical model which is valid for the 
normal daily variations can be accepted for its 
interpretation. The decrease of equivalent height 
and the change of the reflection coefficient of the 
lower ionosphere during a solar flare changes the 
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Figure 5. Monthly means of phase and amplitude of electric 
field strength of GBR during sunrise measured in Berlin in 
June 1962 and December 1962. 

The time difference M is centered to the time of sunrise in 35 km height above 
the middle of I he propagation path. 

interference pattern between groundwave and sky- 
waves. The change of phase, the well-known sudden 
phase anomaly (SPA), generally is positive. The 
change of amplitude can become positive or negative 
depending on frequency, distance, and time of day 
[Volland, 1959]. 

Figure 6 gives an instructive example of three 
successive solar flares on GBR measured in Berlin. 
All three effects are recognized uniquely as solar 
flare effects by other independent, simultaneous, 
solar-terrestrial events and show an increase of 
phase indicating a decrease of equivalent height. 
The amplitude changes otherwise are positive, nega- 
tive, and zero, because the mutual phase of the 
groundwave and the one hop wave becomes slightly 
negative, positive, or zero. 

3. Interpretation of Phase Measurements 

3.1. Ionospheric Model 

The results of phase measurements mentioned 
above can be explained in their essential features 
with the help of VLF propagation theories. These 
theories presume that the region between earth's 
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Figure 6. Amplitude and phase of GBR field strength observed 
in Berlin during three successive solar flares. 

The arrows indicate the beginning of the effects. 



surface* and ionosphere is a plane or curved wave- 
guide with sharp bounded homogeneous walls. 
Electromagnetic energy in such a waveguide can 
propagate only in discrete 4 modes. Solely for short 
distances, the familiar ray theory of wave propaga- 
tion is a suitable approximation. The transition 
between ray theory and mode theory occurs at 
medium distances for VLF waves. Both theories 
give the same results in the transition region 
[Volland, 1961]. 

The diurnal phase lag depending on distance can 
be explained for greater distance's according to Wail 
[1961], Westfall [1961], and Crombie et al., [1958] as 
the difference between daytime phase velocity and 
nighttime phase velocity of the first mode. This is 
caused by the change of equivalent height of reflec- 
tion from about 70 km at daytime to about 90 km at 
nighttime. Thereby the curvature of the earth's 
surface is important, since within a curved wave- 
guide, the effective distance between transmitter and 
receiver changes with the height of the waveguide. 
In this connection, Wait [1963] explained, qualita- 
tively, the phase lag minimum in figure 2 as the 
interference between first and second mode. On the 
other, the frequency dependence of Ar m /p in figure 4 
may be interpreted in terms of frequency dependence 
of the propagation constant of the first mode 
[Blackband, 1962b]. In fact, Wait [1963] explains 
the measured minimum in figure 4 by the exact mode 
theory of curved waveguides. 

In this section the ionosphere is considered as 
plane, but anisotropic and inhomogeneous in vertical 
direction. Concepts like equivalent height, effec- 
tive conductivity, and effective dielectric constant 
will be connected with realistic ionospheric models. 
Since most of the phase measurements have been 
made with the transmitter GBR (Rugby), the 
calculations will be restricted to ionospheric models 
in medium northern geographical latitudes and to 
the frequency 16 kc/s. We expect results which 
are at least typical for VLF waves outside the 
equatorial regions. Detailed calculations of VLF 
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propagation valid for the equatorial zones have been 
made by Crombie [1960, 1961]. 
We assume a Chapman profile, 

N=N m exp{(l/2)(l-exp [(z m -z)/H])} , (16) 

taking into consideration a recombination coefficient 
decreasing exponentially with height and a collision 
number profile according to Nicolet [1959] which 
behaves approximately as 

v=v Q exp [(z —z)/H]; ^0=70 km 

H=8 km (17) 

below 90 km. We neglect the difference between 
v and the effective value of v at low frequencies 
which may differ by a factor 2 [Ginzburg, 1961], 
because already the data of Nicolet are rather un- 
certain, and v itself surely varies with location and 



season. The earth's magnetic field is assumed to be 
vertical with the value of the (angular) gyromag- 
netic frequency at medium northern latitudes 

co ff =8.10 6 s" 1 . 

For this model, and a set of parameters N m and 
z m , the elements of the reflection matrix 



R = 






(18) 



of plane electromagnetic waves depending on angle 
of incidence have been calculated [Volland, 1963]. 
Since the penetration depth of VLF waves at day- 
time is of the order of one wavelength, the ionosphere 
above the D layer has only a small influence on the 
reflection characteristics and can be neglected. At 
nighttime the VLF waves can penetrate into the 
higher ionosphere. The calculated reflection co- 
efficients of the nighttime ionosphere, therefore, 
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Figure 7. Calculated amplitude {fig. 7a) and phase (fig. 7b) of the elements of the reflection matrix for different ionospheric 

models. 
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are more uncertain than the values of the daytime 
ionosphere, especially the phase. 

The results of calculation of the elements of R 
for the parameters JV OT =10 3 cm -3 and z m =75, 85, 95, 
and 105 km are shown in figure 7a (amplitude) 
and figure 7b (phase) . Because of the special choice 
of the direction of the geomagnetic field, j_i?|| = ||.R_L. 
The elements of R are represented as 



R=\R\e"*++\ 



(19) 



In figure 7b, # is plotted versus 0=cos f. 

From the phase <j> an equivalent reflection height 
h' can be defined [Budden, 1961a]: 



,,_ ±f_z#_\ 

h - z <« 2i\dco 8 f/ r=r ; 



(20) 



z vel is an arbitrary reference height of the calculated 
phase. In our example, it is 



Zret=Zm — 20 Kill. 



(21) 



Thus, it is inferred that spherical waves at angle of 
incidence f 0? are reflected at height h/ with a reflec- 
tion coefficient of amplitude \R\ and a phase jump 
of 180°. It is to be noted that h' depends on the 
angle of incidence f. In figure 8, h f calculated from 
the data in figure 7b in conjunction with (20) is 
plotted versus z m for oblique incidence (cos f<0,3) 
and vertical incidence (cos f=l). The slope of h' 
is different at vertical and oblique incidence for the 
same ionospheric model. Qualitatively, the same 
phenomenon has been observed by Bracewell et al. 
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Figure 8. Equivalent height h' versus height of Chapman- 
layer Zm deduced from the ionospheric models in figure 7b 
for vertical incidence (C = l) and oblique incidence (C<0.3). 



[1951] for the transition from night to day. As can 
be seen from figure 8, BracewelPs interpretation of 
the phenomenon [Bracewell and Bain, 1952] assuming 
two different D layers seems to be unnecessary. 

At the oblique incidence (cos f<C0,3), it is seen 
that 

^-^-r« constant. (22) 

dcos f 

This means that h' is constant for cos f<0,3. Now, 
cos f =0,3 is equivalent to a distance of p=450 km 
for A' = 70 km. Therefore, it is possible to operate 
with a constant equivalent height in the whole region 
of medium and large distances. 1 nis fact facilitates 
the introduction of an isotropic and homogeneous 
equivalent of the ionosphere which behaves like the 
inhomogeneous and anisotropic ionosphere for plane 
waxes at oblique incidence and a given frequency. 
This equivalent has an effective conductivity o- and 
an effective dielectric constant e, so that \\\R\\\ and 
|±/«\l| are the same for the inhomogeneous iono- 
sphere and the surrogate at cos J*<0,3. Since the 
phase of the reflection coefficient of a homogeneous 
ionosphere with sharp boundary al oblique incidence 
differs from 180°, the equivalent ionosphere must 
have the virtual height h which is slightly different 
from W [Volland, 1964a]. For \\R U the statement 



\\R±~yC# 



(C<0,3) 



(23) 



is appropriate. The phase term in (23) relates the 
phase of ijJSjl to the height h of ||B||. In figure 7a, 
some calculated reflection coefficients of the equiv- 
alent are drawn as dashed lines which fit to the 
corresponding values of the models with the param- 
eter 2^=105 km. All parameters of the equivalent 
of the models in figure 7 are put together in table 1. 
We see that e can take negative values. 



Table 1. Parameters of equivalent ionospheres which at oblique 
incidence and 16 kc/s have the same reflection characteristics 
as the inhomogeneous and anisotropic ionospheric models (h',h, 
and z m in km) 
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In this manner each arbitrary ionospheric model 
can be substituted by a homogeneous ionosphere 
with a sharp boundary which has the same reflection 
characteristics as the model at oblique incidence and 
a given frequency. 

The influence of the geomagnetic field increases 
with increasing height z m and with increasing cos f . 
At (7=1, even for low values of z m (daytime condi- 
tions), it yields |||flii| = LB±|»|||Bj_|. The reflected 
waves, therefore, have nearly circular polarization. 
With decreasing (7, the value 



InBlI 



hifinl 

decreases. This is the reason for the small influence 
of the geomagnetic field on VLF at oblique incidence. 
Since h and W are constants at oblique incidence, 
the equivalent ionosphere can be applied for the 
calculations at mode theory, too. The effective 
parameters of the equivalent ionosphere A', <r, and e 
strongly depend on frequency [Volland, 1963]. 
Therefore, it is necessary to repeat the calculation 
of R for the same model at other frequencies. 

3.2. Propagation Formula 

The theory of propagation of VLF waves in plane 
or curved waveguides with sharp boundaries has 
been treated extensively. Detailed information may 
be found in the books of Budden [1961b] and Wait 
[1962]. In the VLF range, the transmitter antennas 
mostly behave like vertical electric dipoles. In what 
follows it is assumed that both transmitter and re- 
ceiving antennas are situated on the ground. 

The vertical component of the electric field strength 
of a vertical electric dipole according to a simplified 
form of mode theory is [Volland, 1964a] 



Sl'*coeh*A 1 .e*( 1 " S ""S*')-^ 



(24) 



with 



c 



z Cn + TT A w (eigenvalue of the nth mode) 



^_>-1/2)tt 



A„ 



kh 

TliC n 



1 , yH 12 n 4 Ki C„H— 4- 
e n,C, 






\ € CO/ 

\ e w / 



ja \ /2 (refraction index of the homogeneous 
equivalent ionosphere) 



2 (refraction index of the homogeneous 
earth) 



e, a, h = effective parameters of the ionosphere 
<r e = conductivity of the earth 
w = angular frequency 
6= angular distance between transmitter and 

receiver 
a = earth radius 
^0= undisturbed electric field strength of a ver- 
tical dipole in free space and distance p. 



The mode factors 
B n = 



V: 



jhkp 



S n 



sin 6 



(25) 



take into consideration the curvature of the earth. 
Since S n ~l, B n can be drawn before the sum in (24) 
as a common factor: 



B= 



- /x 

V sin 6 



jhk 



(26) 



For a first approximation, it is sufficient to neglect 
the imaginary part of S n in (24). The phase of 
the nth mode then becomes 



<Pn~ \ 



■sM 



(27) 



In this formula, it is remarkable that the phase 
retardation of the nth mode relative to the un- 
disturbed propagation is compensated by the in- 
fluence of the curved earth on account of 

C 2 n ~hla 

for low order modes. This means that the phase 
velocity of the first mode 



0i= 



'-il^-iM 



(28) 



which is greater than the velocity of light c at 
frequencies below 8 kc/s [Jean et ah, 1960] can be- 
come smaller than c at greater frequencies [Wait 
and Spies, 1961]. A more exhaustive discussion 
of this phenomenon is given by Wait [1962] who 
employs higher order approximations. He indicates 
that the equation is qualitatively correct. 

The difference between daytime and nighttime 
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Figure 9. Polar diagram of E t /2E B for the daytime model depending on distance. 
(a) The left hand diagram is calculated by ray theory, (b) The righl hand diagram is calculated by mode theory. 



phase of the first mode, according 
theory, is 



to the simplified 



knight <Pdi 



L n f *L (A 1 \ | ^ni,ht— A day \ 

P \8k 2 \hl y KUtP ^ / 



(29) 



The influence of the geomagnetic field is placed 
into the values e, a, and y. Now the magnitude of 
the third term of A n is very small compared with 

the other terms on account of (7f/4<< 1 JU1( 1 ma . v 
be neglected. The consequence is that in the ray 
theory approximation the terms with \\R± and 
±E\\ also may be neglected. The corresponding 
formula of ray theory, therefore, is 

E 2 ^2E q b(w+2 ± sin 3 f,. Q^J 



fir 1 lii#i'lv* ( '- r n ) -.'»* 



(30) 



where the Sommerf eld function Wis replaced by unity, 
^n=— , r n =-yJp 2 +(2nh / ) 2 ) A' = equivalent height, 



sin 



\\\H\\ | = amplitude of the reflection coefficient of the 
inhomogeneous and anisotropic ionosphere, and, 
finally, R e = reflection coefficient of the earth for 
vertically polarized waves. 

The curvature of the earth has been taken approxi- 
mately into consideration by the factor B (see (26)). 
If we assume the earth's surface to be ideally con- 
ductive (A en =0; B e = l), the resultant error pri- 
marily influences the amplitude of E z . In the 
following calculations we shall make this assumption. 
The neglect of the finite conductivity of the earth, 
however, is allowed for only up to frequencies of 
about 20 kc/s. 



3.3. Interpretation of the Diurnal Phase 
Measurements 

With help of the propagation formulae (24) and 
(30) together with the ionospheric parameters given 
in table 1, the field strength E z has been calculated 
depending on distance p. As the daytime model, 
the model with the parameters 7V W =1000 cm" 3 and 
2 m =85 km (h' = lQ km) was selected. The reason 
for the choice of (his model may be seen in section 3.5. 
For the nighttime, the model with the parameters 
iV m =1000 cm" 3 and 2^=105 km has been used. For 
convenience, an equivalent height of /*/ = 90 km 
instead of 88 km, according to table 1, has been 
introduced. 

The calculations have been made by (30) up to 
p=1500 km and by (24) above p=1000 km. The 
agreement of the results in the range between 1000 
and 1500 km served as control. The calculated 
difference — A<p m is shown in figure 2 as a continuous 
full line. Comparing the calculated curve with the 
experimental data, we see that an agreement is pos- 
sible only if we take two jumps of 360° at p=600 km 
and p=2600 km. This discrepancy can be explained, 
if we plot E Z /2E B in a polar diagram. This is 
shown in figure 9 for daytime propagation and in 
figure 10 for nighttime propagation. Figure 9a 
originates from the ray theory (30). The straight 
line 0^4 corresponds to the groundwave (assuming 
1^=1). The thin straight lines — drawn in figure 9a 
only for the case p=T500 km — correspond to the 
multihop skywaves. The portions of the skvwaves 
increase with increasing distance. The vector drawn 
from the point to an end point of the curve gives 
amplitude and phase of EJ2E B. Figure 9b is 
originated by mode theory (24). Here the straight 
lines correspond to the mode terms. The modes of 
higher order decrease with increasing distance. 
Figure 10 must be read in the same manner. 
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Figure 10. Polar diagram of E Z /2E B for the nighttime model depending on distance. 
(a) The left hand diagram is calculated by ray theory, (b) The right hand diagram is calculated by mode theory. 



From figure 9a and 10a, we see that between 500 
and 1000 km the origin will be encircled by the 
endpoint of the vector in a clockwise sense at day- 
time and in a counterclockwise sense at nighttime. 
The nighttime phase, therefore, receives a gain of 
360° or one wavelength in respect to the daytime 
phase. We would obtain such a result if the meas- 
urements were made by a synchronized oscillator in 
an aircraft flying off the transmitter at daytime and 
nighttime, respectively. Yet, there exist only few 
measurements of this kind [Reder et ah, 1961]. All 
measurements plotted in figure 2 are made at one 
location. By this fact the phenomenon is hidden. 

Figure 11 shows two polar diagrams for two con- 
stant distances, p=600 km and p=650 km. Here it 
was assumed that the equivalent height h' and the 
corresponding ionospheric parameters in table 1 
changed continuously from the nighttime value to 
the daytime value. We would expect such picture 
if we would measure continuously, phase and ampli- 
tude at one location. As can be seen from figure 11, 
the origin is circulated in a left hand sense at the 
smaller distance and in a right hand manner at the 
larger distance. There must exist a distance between 
both these distances where the vector runs through 
the origin. The phase becomes indefinite. In this 
range, the measured diurnal phase lag jumps by 360° 
from negative to positive values. This is indicated 
in figure 2 by the vertical dashed line at 600 km. 
Receiver stations situated round this critical distance 
may fall from the positive into the negative phase 
lag region, or vice versa, if the reflection character- 



istics of the nighttime ionosphere changes from one 
night to the next. The phase difference must then 
differ rather exactly by 360°. This is the qualitative 
interpretation of the phase anomaly observed in 
Lindau (see fig. 3). 

In our model, a second phase gain of 360° of the 
nighttime phase occurs at 2600 km. Here the sum 
of the second and higher order modes is larger than 
the first order mode and is phase shifted by 180°. 
The magnitude of the minimum of - A<p w in this 
region depends very sensitively on the minimum 
difference vector of the modes at nighttime. Un- 
fortunately, the measurements in these distances are 
insufficient to decide whether the second phase gain 
is real or not. In any case, phase measurements in 
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Figure 11. Polar diagram of E Z /2E B depending on equiva- 
lent height h' for two different distances p = 600 km and 
p = 650 km. 
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these two regions give very sensitive indications of the 
diurnal and interdiurnal change of the lower iono- 
sphere. 

The agreement between measurements and calcu- 
lation in figure 2 can be only of qualitative kind be- 
cause the measurements are made at different times, 
different regions and different frequencies. The 
dash-dotted line in figure 2 is the calculated phase 
lag of 1 he first mode, according to (29). We see that 
up to 5000 km the second mode may not be neglected 
for the nighttime propagation. 

The value - Ar m /p calculated from (29) for the 
first mode is shown in figure 3 as a dashed line. Here 
a nighttime height, A'night.^88 km, has been chosen. 
The theoretical curve fits well with the measured 
values. No attempt has been made to explain the 
minimum at 16 kc/s, since the number of measured 
frequencies in the region is very small. 

3.4. Interpretation of the Sunrise Effects 

The model previously used is valid only if it is able, 
also, to explain the transition from night to daytime 
conditions. At shorter distances, when the one hop 
sky wave exists primarily besides the ground wave, it 
seems to be appropriate to approximate the transition 
by lowering the whole ionosphere continuously. In 
figure 12a, calculated values of EJ2E B, according to 
ray theory, are plotted in the usual manner for the 
distance p=980 km and different heights h'. Ampli- 
tude and phase from figure 12a are drawn separately 
in figure 12b as full lines versus h'. Here the phase 
term of B has been added to the phase. The dashed 
lines in figure 121) are amplitude and phase of GBK 
measured in Berlin on December 12, 1962. This 
sunrise effect represents well the monthly means in 
figure 5. The scale of W of the theoretical curve and 
the time scale of the measured values are related in 
such a manner that the two curves fit best. The 
sunrise effect begins at a zenith distance of x=99° 
at the middle of the propagation path and ends at 
X = 93.5°. The agreement between calculation and 
measurements is very good. It is worth mentioning 
that the scaled h f must not be proportional to the 
registration time . 
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Figure 12a. Polar diagram of E Z /2E B depending on equiva- 
lent height h' for p=980 km. 



The sunrise effect in summer (see fig. 5) cannot be 
explained as well as the winter effect. But we get a 
qualitative agreement between calculation and 
measurement, when we assume that in summer the 
nighttime ionosphere ends at h' = 85 km. Then the 
deep minimum of amplitude and the turning point of 
phase caused by the interference between one hop 
and two hop wave remain, but the amount of 
amplitude is too small in respect to the measured 
values. If we suppose a larger amplitude of sky- 
waves in summer, we even get quantitative agree- 
ment. 

A likewise qualitative agreement exists between 
measured and calculated sunrise effect at distances 
of 2150 km. In figure 13, a diurnal phase change of 
GBR measured in Malta (p=2150 km) by Black- 
band [1961] has been plotted together with calculated 
data according to our model for the transition from 
night to day. The course of phase is the result of 
interference between first and second mode. We 
cannot expect better agreement , since the validity of 
our transition model begins to cease at these dis- 
tances. 

3.5. Interpretation of Daytime Phase variations 
and Solar Flare Effects 

The normal variation of equivalent height at day- 
time is caused by the variation of the electron 
density profile of the lower ionosphere. During a 
solar flare, the illuminated part of the lower iono- 
sphere becomes additionally ionized, and therefore 
the equivalent height is lowered for a short time. 
If one assumes that the region of the lower ionosphere 
where the reflection of the VLF waves takes place 
behaves approximately like a monochromatic Chap- 
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Figure 12b. E Z /2 E from figure 12a, amplitude and phase 
separated, versus equivalent height h' (full line). 

The dashed lines are amplitude and phase of GBR field strength during sunrise 
versus time observed in Berlin on December 12, 1962. The arrows indicate 
time of zenith distance x at the middle of the propagation path. 
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Figure 13. GBR phase during sunrise versus time observed 
in Malta on May 20, 1960, according to Blackband [1961] 
{full line). 

The dashed line is the calculated phase versus equivalent height h'. 

man layer, then the daytime variation of the electron 
profile is caused by the variation of zenith distance 
X of the sun. On the other hand, during a solar 
flare, the additional solar emission at nearly constant 
zenith distance can be expected to increase the 
strength of the layer without changing its shape. 

Recently, however, this conception has been 
questioned by Friedman [1960], who indicated that a 
new layer caused by X-rays is built up during a 
solar flare. In the following, only the first conception 
of one layer will be used. For the second conception 
of two layers [see Volland, 1964b]. 

Ionospheric reflection coefficients V R\ { have been 
calculated for 16 kc/s and a number of ionospheric 
models with electron profiles like (16) and collision 
number profiles like (17) [Volland, 1963]. The 
chosen parameters were z m = lb and 85 km and N m 
between 100 and 10 5 cm -3 . From this calculation it 
follows that the amplitude \\\R\\\ does not change 
very markedly (see the values of | ^R^ | in fig. 7a for 
z m = 75 and 85 km). On the other hand, the phase 
and therefore h' show essential variations. 

A sharply bounded homogeneous surrogate iono- 
sphere with parameter n\=\—j has nearly the same 
value |||i?u| as the inhomogeneous models at oblique 
incidence. The imaginary part of ri\ is 



6 CO 



(31) 



Moreover, we state the well-known relation for very 
low frequencies 



W 



N 
2.81. 10" 8 — mho/m 



(32) 



(e and m are charge and mass of electron) . 

We now fix that height z a within the inhomogene- 
ous ionosphere where (31) and (32) yield: 



N{z a ) 



e o; 



v(z a ) 2.81 



• 10 8 =31.5m- 3 s 



(33) 
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Figure 14. Relation between equivalent height h r and the 
height z a of the point N = 31.5i> of a Chapman-layer in the 
height z m at frequency f = 16 kc/s. 

Points and circles are calculated from the reflection coefficients \\R\i of different 
models. The dashed straight line is 

2 m — Z a Zm — h' 
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and calculate the difference z m —z a . If we plot this 
value versus the difference z m —h f of the same model, 
we get figure 14. The points gained from the dif- 
ferent models nearly lie on the straight line 



-h' 



H 



H* 



(34) 



with H=S km (scale height), fl r *=6.35 km (equiva- 
lent scale height), and z* l = z m — 5.7 km (equivalent 
characteristic height). If (z m —z a )/H is replaced by 

(3*-A')/ff*in 



N(z a )__N m 

v(z a ) V m 



exp 



{it 1 



^\Zm % a) 



— exp 



H 

(Zm 



•Za 



H 



]}■ 



=31.5 (35) 



it is seen that each measured value of h' corresponds 
to a height z a where A^=31.5^. Since v is supposed 
to be known, we are able to determine N in the 
height z a from the measured value of h f . For 16 
kc/s and oblique incidence, the lower ionosphere in 
daytime, therefore, behaves like a sharply bounded 
homogeneous ionosphere with refractive index n 2 
= l—j. Its height h' is in the point Af=31.5^ of a 
transformed Chapman layer with the parameters 
N m , z* and H*. 
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In view of the relation 

z m =z mn —H In cos x 



(36) 



where z m0 is the characteristic height of the Chapman 
layer at zero zenith distance, it follows from (35) 
that 



In 



1 



h'-s* 



COS X 



H* 



+ln/: 



-2 In 31.5 ^+2 

N m0 



h '-fA \. 



(37) 



This equation makes possible the determination of 



z m0 and 






if the daytime variation of h' is known. 



Figure 15 shows values of h' determined from the 
field strength of GBR measured in Berlin between 
April and June 1959 [Volland, 1959]. The dashed 
line is calculated from (37) with help of the parame- 
ters: z m0 =S5 km, N mQ =l0 9 m- 3 =1000 cm -3 , v m0 = 
1.5-10 6 s~ l in 85 km height (see (17)). At the 
reflection point of the one hop wave at local noon in 
May, it is cos x = 0-82. From this it follows that 
m =86.5 km and iV m =iV w o=100() cm -3 which is in 
good agreement with the daytime model hitherto 
used. 

During a solar flare, the parameter N m is changing 
at constant z m . The retardation of N in respect to 
the solar emission is governed by the equation 



(lN •ii*- 



)-a 7 V m 



(*£0) 



(33) 



for m=l or 2. Here g is the normal ion production 
rate prior to the onset which is considered to remain 
constant during the flare; 5 is the maximum relative 
increase of ion production rate; e~ yt is the ideal time 
course of a typical solar flare; and y is the decay 
constant. In (38), m = 2 corresponds to the effective 
recombination in which case a is a recombination 
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Figure 15. Equivalent height h' determined from field strength 
measurements of GBR in Berlin between April and June, 
1959, versus time {points) and theoretical curve (dashed 
line) . 



coefficient. On the other hand, m=l corresponds to 
effective attachment. In which case, a is an effective 
attachment coefficient. 

Neglecting vertical motion of electrons and pre- 
suming values of y, 8, a, and m, the momentary 
relative increase of electron density 



AN(z, t) 



(39) 



can be calculated from (38) for each height z, de- 
pending on time At from the beginning of the effect. 
The height z a can be determined from these electron 
profiles and from (33). Likewise the difference 
Ah'=h'(t)—h' () can be calculated from (34). Other- 
wise Ah (t) can be determined from measurements 
during a solar flare. Figure 16 shows as an example 
Ah f versus At. The points arc measured values of a 
solar flare effect. The full line is evaluated with 
help of the parameters: 



5 = 85 
7=7.2- 10-V 
a=2.2.10"V 
m = 1 . 

Here it was presumed 



1 in 70 km height 



a oc e 



-z/H 



(40) 



(41) 



The combination of parameters (40) is unique, since 
the curve must fit both height and time of maximum 
and the decay of the measured values with time. 

Similar calculations have been done for 15 additional 
selected solar flare effects [Volland, 1964b]. The 
selection was taken in such a way that only solar 
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Figure 16. Course of departure of equivalent height from the 
normal value Ah' during a solar flare determined from field 
strength measurements of GBR in Berlin (points) and 
theoretical curve (full line) versus time after beginning of the 
solar flare. 
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flares were used which had the proposed course of 
emission according to (38). To decide this question , 
the 3.2 cm solar radio emission which is supposed 
to behave very similarly, has been used. Twelve of 
the fifteen effects could be interpreted quantitatively 
by (38) only with the parameter m= 1. The medium 
value of the attachment coefficient is 

a=3.4.10-V 1 in 70 km height 

in agreement with results of HF-absorption meas- 
urements during SID's [Taubenheim, 1962]. The 
rest were very strong effects with maximum 
height differences Ah' of more than 9 km. The 
effective attachment law therefore seems to be 
valid in this height. 
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